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Abstract 

We characterize quantum oscillations in the magnetic susceptibility of a quantum 
critical non-Fermi liquid. The computation is performed in a strongly interacting 
regime using the nonperturbative holographic correspondence. The temperature de- 
pendence of the amplitude of the oscillations is shown to depend on a critical expo- 
nent V. For general v the temperature scaling is distinct from the textbook Lifshitz- 
Kosevich formula. At the 'marginal' value = i, the Lifshitz-Kosevich formula is 
recovered despite strong interactions. As a by-product of our analysis we present a 
formalism for computing the amplitude of quantum oscillations for general fermionic 
theories very efHciently. 



1 Results and background 



A central question in the theoretical characterisation of non-Fermi liquids is the fate of the 
Fermi surface. For instance the 'strange metal', and perhaps quantum critical, regions of 
the cuprate or heavy fermion phase diagrams separate phases with very distinct energy- 
momentum distributions of fermions. This is seen in many experimental probes, a recent 
discussion with strong overlap with the concerns of the present paper is [1]. 

Quantum oscillations are a robust feature of systems with a Fermi surface [2]. The 
recent and ongoing experimental observation of quantum oscillations in the copper oxide 
high temperature superconductors [3llll[5l|6l[71[8l|9] is reinvigorating theoretical approaches 
to the subject (e.g. [IHIIII])- Present measurements are, perhaps surprisingly, consistent 
with textbook results for quantum oscillations in Fermi liquids. However, as new regimes are 
investigated, in these and other quintessentially non- Fermi liquid materials, it will be crucial 
to have theoretical templates available for comparison. For instance, the exciting recent 
results of [12] show that the effective quasiparticle mass, as read off by fitting quantum 
oscillations to the Fermi liquid formula, appears to diverge as one approaches a metal- 
insulator quantum phase transition. A similar divergence is observed in heavy fermion 
compounds [13j. Here one encounters a theoretical hurdle; the most interesting regimes 
are often strongly coupled, and perturbative quantum field theory treatments may not 
fully capture the physics of interest. Recent work investigating the effect of interactions 
on quantum oscillations includes [TU [151 [El (El- In this paper we will use the inherently 
non-perturbative 'holographic correspondence' (see e.g. [HI [19] for relevant introductions) 
to give a controlled computation of quantum oscillations in the magnetic susceptibility 
of a strongly interacting quantum critical non-Fermi liquid. We will however highlight 
similarities with the approach in [T^ . 

The main result of this paper will be the following expression for the leading period de 
Haas - van Alphen magnetic oscillations in a class of 2+1 dimensional theories that exhibit 
an emergent quantum criticality at low energies: 
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where x is the magnetic susceptibility, e is the charge of a fermionic operator, A is the area of 
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the sample, T the temperature, c the speed of hght, kp the Fermi momentum, B the apphed 
magnetic field, /x the chemical potential and u a critical exponent. Our computations are 
in the clean limit, with no disorder. The most important of these parameters, for our 
purposes, is the critical exponent u, which satisfies 0<z/<|. Atz/=^we will find 
Fn{^) = '2n'^h{n + |), where h is a. dimensionless constant defined below. The sum in 
expression ([T]) then gives 
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This is essentially the textbook Lifshitz-Kosevich result [21 [20], as we discuss in more 
detail below. Our theories will be in 2+1 dimensions, although many results can likely 
be generalized to 3+1 dimensions. When z/ < ^ the functions -Fn(^), given below, are 
considerably more complicated. The point we wish to emphasize, however, is that at larger 
temperatures T > the decay of the amplitude as a function of T is not of the simple 
exponential form predicted by the Lifshitz-Kosevich formula, but rather 

Xosc. ~ e"^'" . (3) 

This is what we will mean by a generalized Lifshitz-Kosevich scaling. If we write this scaling 
as a temperature dependent effective quasiparticle mass in the usual Lifshitz-Kosevich 
formula, then 

~ — U ' (4) 
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which is divergent when T ^ /i and u < ^. This is perhaps interesting in the light of the 
observations in [121 US] . 

While the large temperature scaling (|3| is the most universal feature of our results, we 
can also plot the full amplitude ([T]) as a function of temperature for given values of the 
parameters. We will introduce the various free parameters of the model below. Typical 
results are shown in figure [Tj The most interesting observation is that for a given value of 
the critical exponent v < \ there is a range of possible behaviors at low temperature. While 
the curves can saturate, mimicking the usual Lifshitz-Kosevich behavior, it is also possible 
for the curve to reach zero temperature with a finite negative gradient or alternatively to 
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exhibit a maximum before reaching zero temperature with a positive gradient. A maximum 
was reported experimentally in [21]. In [21] it was further noted that an improved fit to 
the data could be achieved by modifying the Kosevich-Lifshitz formula. 
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Figure 1: Typical dependences of the amplitude of quantum oscillations on temperature. For 
illustration = 5, = 1, a = 1. Angles of h from top to bottom: ip = {—(fo, —0.2(fo, 0.51(/?o, 
(fo} where the maximum value ipo = tt{^ ~ The magnitude of h has been scaled to make the 
large temperature behavior coincide: h = {0.34, 0.39, 0.58, 1}. 

The theories for which ([T]) will be shown to hold are described using the 'holographic 
correspondence'. We will not review the methodology in detail, introductions written 
for the condensed matter community can be found in [181 HH] 5 but rather summarise the 
physical properties of the theories in question. 

The holographic correspondence allows a class of strongly interacting quantum field 
theories to be studied in a limit in which there are a large number of degrees of freedom 
per site. Unlike more traditional vector 'large limits, the theories do not become weakly 
interacting in this limit, and might therefore be expected to capture aspects of interesting 
experimental systems that would otherwise elude theoretical control. 

It was shown in [22], following earlier work in [23l [Ml ES], that the fermion spectral 
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densities in these theories exhibit a broad peak with a zero temperature dispersion relation 
at /c ~ fci? of the form 

— + he'^u^"" = k-kF, (5) 
vf 

where {vF,h,9,h',kF} are real constants. For < ^ the non-analytic term dominates 
at low frequencies, leading to non-Fermi liquid behaviour. The case = ^ leads to the 
dispersion relation ^ + he'^^ulogu = k — kF, which is precisely that of a marginal Fermi 
liquid [26j. For all z/ < |, the peak in the spectral density does not correspond to a stable 
quasiparticle excitation. This is because the width of the peak is always comparable to its 
height. Viewed as a pole in the spectral density in complex frequency space, its residue 
goes to zero as the pole hits the real axis at k = kF [22] • In principle we could also study 
u > ^, but here the linear term in ^ dominates at low energies and a more conventional 
behaviour is expected. See however |22] for some curious properties of these cases. 

Given that ^ does not describe a weakly interacting (stable) quasiparticle, one can 
anticipate that the contribution of the fermions to thermodynamic and transport quantities 
will not be simply that of a free fermion with dispersion The correct way of computing 
in these systems was developed in |27j , with the more mathematical aspects treated in [28] . 
The essential step is to consider ^ as the singular locus of the fermion spectral density 
ImG/j(w, fc). It is easy to see that ([s]) has two types of singularities, a pole and then a 
branch cut emanating from u = 0. While the pole describes the naive 'quasiparticle', both 
the pole and the branch cut will give contributions to e.g. thermodynamic quantities. 

This paper will be concerned with small but finite temperatures. At finite temperature, 
the branch cut of (|5| is resolved into closely spaced poles. For T,cu <^ fi one obtains 
that the poles of lmGji{u, k) are given by solutions to 



J-K(A;)) = 0, (6) 



where 
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See e.g. figure (3) of [27]. The dimensionless constant a is related to the normalisation 
of the current-current correlator [22]- While complicated, this formula is largely fixed by 
an emergent SL{2,'R) (or possibly even Virasoro) symmetry at energies u <^ fi, suggesting 



perhaps validity beyond the specific holographic theories considered in [22]. This emergent 
IR scaling symmetry is the quantum criticality referred to in the title of this paper. The 
only dimensionful scales in the theory are the chemical potential fi, magnetic field B, Fermi 
momentum kp and temperature T. In ([t]) we have assumed that u < ^ so that the linear 
in u term in ^ can be dropped at low energies. 

All of the poles given by ^ contribute to quantities of interest, even those that are a 
long way away from the real frequency axis. The key result of [271 12H] was to express the 
contribution of the fermions to the free energy as a sum of contributions from these poles. 
The formula is 

„R4r /, ^-.(O^yy 
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Anticipating our interest in magnetic fields, we have given the free energy as a sum over 
Landau levels rather than momenta. The first term in ([s]) is the degeneracy of the Landau 
levels. The frequencies c<;^,(£) are obtained from in ([6]) and ([T]) by the replacement 

^2 _^ This replacement is precise in the limit ^ <^ kp that we will be interested 

in. The formula ([s]) is not as exotic as it may appear; for instance, the free energy of 
a damped harmonic oscillator can be computed using essentially the same formula, with 
o;^ again given by the poles of the retarded Green's function [271 EH]- The appearance of 



\T(ix + is a generalisation of the Fermi-Dirac distribution to complex energies. If x is 
real then \r{ix + = Trsechvra;, recovering the standard expression. 

Our objective is to perform the sum ([s]) given ([t]) to obtain the magnetic susceptibility 
for general T ~ <^ /i. The result for the leading oscillatory part of the susceptibility is 
stated in ([T]). 



2 The computation 

Our starting point is the formula for the fermionic contribution to the free energy, given in 
^ in terms of the poles ^ of the fermion retarded Green's function. It will be useful to 
consider the dimensionless quantity 

^ ^ = Re J] J] 21ogr + I) , (9) 
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where we set 



In the formula ([T]) defining the poles we will furthermore set 

n = = Aifsm (p + I cos ip) , ( 11 ) 

so that {h, h, (p} are now dimensionless. While in principle these parameters are determined 
by data in the UV by solving some ordinary differential equations numerically [22], we will 
simply treat them as order one quantities, as we are more interested in parametrising 
possible low energy physics. There is a restriction on ip that ensures that the poles are in 
the lower half frequency plane: — vr(| — u) < p < — u) . Notice that the imaginary part 
of h is always positive. 

Using all these expressions we can rewrite the sum over Co'*(£) as a contour integral. 
Noticing that J-'(x) does not have poles, just zeroes in the right half plane (corresponding 
to the poles of Gr^u, k) in the lower half plane) we can write 
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The contour was chosen such that it leaves the poles of to the right and the branch 



cut of logF (x + I) to the left. Implicitly we are also taking the contour to include a large 
semicircle in the right half plane. We will not need to evaluate the contribution from the 
semicircle explicitly, at a later step we will exhange the current sum over poles inside the 
contour for a sum of poles outside the countour (i.e. in the left hand plane). 



We would like to integrate ( 12 ) by parts, but this is complicated by the presence of the 
branch cuts from the logarithmic term. However, the derivative of Cl with respect to the 
magnetic field can be safely integrated by parts to give 

M.|^.Relvr'""..Mi%3£^. ,13, 



We will be interested in considering the periodic behavior in of this expression. Therefore, 
it is of use to Fourier transform the Landau level variable ^. We will perform a Poisson 
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resummation to rewrite (13). The formula we use is 
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It is straightforward to apply this formula to (13), with the Landau levels going over 



0,1,2,.... We obtain 
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where 

G(a;,5,A;) 
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in which we used the explicit form of (7), changed variables to u 

r (I + z/ — -ia — x) 
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and set 
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Equations (15) and (16) appear to involve formidable sums and integrals. However, 



we can now neatly separate out the oscillating and non-oscillating parts of this expression. 
We will deform the contour in such a way that the integral follows a steepest descent path 
of the exponential term. The reason this helps is that the resulting integral is manifestly 
non-oscillating in 1/B. 



We therefore deform the integral in (16) by m — )■ e*!''! -^u. It is crucial to realize here 



that the contour needs to be rotated in opposite directions in the complex plane, depending 
on the sign of k, to guarantee convergence. The only possible obstructions to this contour 
rotation are either a contribution at infinity, which is absent in our case as the integrands 
decay exponentially if the paths are rotated in the correct direction, or if a pole is crossed 



as the contour is deformed. The expression (16) makes manifest that there is such a pole 
at u = 1 + 7rh (j^j Su{x). In the limit of physical interest, T//i — )• 0, this pole is slightly 
off the real axis, for < u < ^, where our formulae are valid. 

The exact position of the pole depends on the phase of h but it is always slightly above 



the real axis (this can easily be checked for the allowed range of values of h and x G — | 
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The upshot is that for negative k we can rotate the contour and get 

G(^,B,-\k\)^^^e-'i in y — ^. (18) 

This contribution is strictly non-oscillating^ in Deforming the contour for positive k 
we pick up a contribution from the pole. Calculating the appropriate residue yields 

G{x,B,\k\) = G,on-osc.{x,B,\k\)+Gosc{x,B,\k\) (19) 

= G.„(x,5, 1^1) + ^ (i 



The first term is non-oscillating and is the same as (18) with various factors of e*'^/^ — )■ 



e ^^/^. We are therefore left with the following oscillating contribution 

Gosc.(a;, k) = Q{k)G,,,Xx, B,\k\). (20) 

where Q{k) = 1 ioi k > and Q{k) =0 for A; < . The k = term is also non-oscillating 
and does not concern us. We have thus performed the first of our integrals, insofar as 
obtaining the oscillating term is concerned. 



The next integral to address is the x integral in (15). We will convert this integral into 
a sum over residues that are outside the original region of integration. That is, to the left of 
the imaginary axis. Doing this allows us to represent the integral as a sum of the residues 
of the poles of tJ-xtt- These are located at — n with = 0, 1, 2, 3, . . . and have minus 



unit residue. Combining this operation with the result (20), our expression (15) becomes 



M_ = Re||f (1 + 'rh{lfsA-h - „))%-Si'('«'-(i)''-^(-^»>)' . (21) 

k=l n=0 ^ ^ 

It is clear at this point that we have obtained sums over terms that both oscillate and 
decay in 1/B. We can now take the physical T/yU — > limit keeping only leading terms 
determining the oscillations and exponential decay. The result is 

= ^f:sin^f;e--^*(^)^^--^^^(-i-'^) . (22) 

k=l n=0 



more formal way to state this fact is that upon Fourier transforming in -g, the transform has poles 
only at imaginary frequencies. 
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This last formula is essentially the result. To compute the magnetic susceptibility x we 
have to reinsert the factors that relate to f2 in nm. Thus 



d'^Vl eAT ~ eBATdM , , 

^ = -fim= — . 23 

oB^ TIC 2nc oB 

For situations of physical interest we have — kp and therefore the leading result comes 



from the second term by acting with the derivative on the sine in (22). Focusing on the 
leading period, the k = 1 term, this gives our main result, that we already quoted in 
equation Q, with 

Fn{iy) = 2^^1111 hS^ {-l-^) ■ (24) 
We also already noted in the introduction that the case u = ^ is special. This is because 



the ratio of gamma functions in (17) simplifies in this case to give -Fn(^) = 27[^h{n + ^). The 
sum over n can then be done explicitly, to yield a result of the standard Lifshitz-Kosevich 
form (g. 

In general, we cannot perform the sum over n in closed form. However, it is simple 
to check numerically that for all allowed values of the parameters, Fn{v) is positive and 
monotonically increasing in n. Therefore at the high temperatures of primary interest we 
can keep only the first term in the sum in ( [22| or ([T]) given by n = 0. This observation 
also implies that the k = 1 term kept in ([T]) has an exponentially larger amplitude than the 
other terms in this regime. Thus we obtain, for general z/ < ^, the non-Lifshitz-Kosevich 
scaling that we quoted in 

3 A general formula for quantum oscillations 

We will now rederive the result ([T| via a slick argument. The argument is quite general 
and we anticipate future applications. The method used is a generalisation of that in |280 
and we will be brief in presentation. 

The statement is that for any fermionic system satisying assumptions to be given shortly 
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^Frederik Denef has been privately advocating this type of generahsation for some months. 
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where the £*(n) are defined as the solutions to 

^(a;„,4H) = 0. (26) 

Here J^(w, i) = defines the singular locus of the fermion retarded Green's function in a 
magnetic field, Gr^oj,^). The fermionic Matsubara frequencies are Un = 2niT {n+ We 
assume for simplicity that there is a unique ii,{n), but it is simple to relax this assumption. 
It is clear that using n7« with fc^ = solving for £*(n) as in (26) and plugging into (25) 



immediately reproduces our previous result (21). 



The class of theories to which the formula ( 25 ) will most directly apply are those where 
the fermionic partition function can be expressed as the determinant of an operator O in 
a thermal Euclidean space. This certainly applies to free theories and to theories with 
classical holographic duals. In the latter case the determinant is in one extra curved 
spacetime dimension, but this does not make a difference to the argument. We assume 
that in a background magnetic field, the eigenvalues of the operator can be labelled by 
the quantum numbers Un and i as well as any others. The type of reasoning in [28] is 
quickly seen to imply that we must have, up to UV contributions that can be dealt with 
systematically but which will not contribute to oscillations, 

Q = -TtrlogO = Re ^ ^ log(£ - 4(n)) . (27) 

The logic that leads to this expression is to separate the eigenvalues of O according to a;„ 
and i. The contribution from positive and negative u;„ to the determinant are complex 
conjugates of each other [271 EH] so we concentrate on the positive Matsubara frequencies. 
For a fixed a;„ we can deform the operator by letting + 7 and then match the zeros of 
the determinant of (9^,7 as a function of 7. Zeros arise whenever has a zero mode. This 
in turn occurs whenever the Euclidean Green's function has a pole at a; = which we 
define to occur at £ + 7 = li,{n). The retarded Green's function is the analytic continuation 
of the Euclidean Green's function from the upper half frequency plane, thus connecting 



with our definition of Tiyi^ £) appearing in (26). Writing det On,-y ~ Ylei^ + 7 — £*{n)) and 



setting 7 = gives (27). 



Poisson resumming (27) using (14) and picking out the oscillatory part of the Fourier 
transform by rotating the contour in different directions for negative and positive k, in a 
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similar way to how we did previously, then directly leads to (25). Only the rotation at 
positive k leads to a singularity contribution giving the oscillating term. 

We now see that the formula ( [25| reproduces known expressions for free fermions. The 
non relativistic, spinless electron (the effect of spin is simply to multiply the answer by two 
in the limit ^ <t: kp) has 

Be 

^non-rel. (w, i) = i - jJ. - U . (28) 

mc 



It is trivial to solve for £^,{n) defined via (26). Plugging into (25), differentiating twice and 
performing the geometric series sum over n as previously leads to 

_ 27rAT/x^m^c ^ fccos(27rfc^) 

B^e sinh (27r^k^) ' ^ ^ 

k=l V Be / 

This is literally the Lifshitz-Kosevich formula |2j in 2+1 dimensions, which we have derived 



rather painlessly. The fact that J-' is linear in £ in (28) makes the steps leading to (27) 
trivial in this case, there is no rewriting involved. 

We can treat the spinless relativistic fermion similarly. In this case 

J'rci. {to, i) = m^c^ + 2Beic - {fx + . (30) 



It is again immediate to solve for i. Use of (25), the limit T <^ /x, differentiation and 
summing a geometric series gives 



Xo 



( cfe^ 

1^^ ^,^X.(^2V.TJ^\ ^^^^ 



253e sinh(7r2fc^) 

k=\ V Bee) 



We used the relation k?p(? = m?c^ — fi'^. In the massless limit (or /i much larger than mc^) 
this expression recovers our result (j2]) for the 'marginal' non-Fermi liquid at = ^ if we 
choose h = 1. 



The expression (25) is essentially the same as a general expression appearing in [H]. In 
|14j the effects of interactions are incorporated into a renormalised self energy for quasi- 
particles whose one loop contribution to the susceptibility is then computed. This is a 
controlled approximation if there are well defined quasiparticles so that higher order cor- 
rections that cannot be absorbed into the self energy are negligible. In the holographic 
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theories studied here the self energy due to strong interactions is captured by the propa- 
gation of the fermions on a nontrivial background spacetime, leading to the singular locus 
([7]). Interactions between these fermions are suppressed by the 'large N' limit in which the 
holographic computations are performed. Therefore, holography provides a controlled set- 
ting in which the self energy can be strongly renormalised to the extent that there are not 
well defined quasiparticles and yet quantities such as the susceptibility can be computed 



with a determinant formula like (25). 



4 Magnitude of oscillations and the Fermi surface 

We need to check that our result could in principle be measured. For that purpose we com- 
pare the order of magnitude of the amplitude of the oscillating part to the non-oscillating 
part. We will pursue this calculation at low temperatures, where the oscillating signal is 
strongest. In this limit, we will see that the oscillating susceptibility strongly dominates 
over the non-oscillating part in the regime of interest ^ <^ fcl^ for | < z/ < |. This domi- 
nance is, of course, not a strict requirement for experimental detection. We first estimate 
the oscillating magnetization. At low temperatures all terms in the sum in ([T]) are impor- 
tant. In fact, the infinite tail of this sum dominates. Therefore, we can replace S^{—^ — n) 
with its n — )■ cxD limit, Sy{—\ — n) — > v?^ . Because the quantity appearing in the sum 
is ^n, we can replace the sum in n with an integral at leading order in ^. Therefore the 
magnitude of ([T]) becomes 



— ^^./'"^^^^ ~^x(^J • (32) 



It is interesting to rederive this last result from a different perspective that makes 
transparent the role of a Fermi surface. At low temperatures the susceptibility is most 
naturally written as a sum over ^, without Poisson ressumming. We can start from the 



expression (13) for M and calculate x by use of (23). As before, we can change the x 
integral to a sum over poles labeled by n. Once again, at zero temperature the tail of this 
sum dominates and we can substitute ^„ f dn and Su{—^ — n) — )■ n^^ . The resulting 
integral can be performed analytically to leave a sum over ^ that is similar to the expressions 
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obtained in [27j. This sum has a nonanalyticity at i = Expanding the susceptibihty 
at small ^ using a generalized version of the Euler-Maclaurin formula |30j, the sum in i 
becomes an integral plus contributions at the edges. The edge near the Fermi surface is 
responsible for the leading effect we are interested in. Explicitly 

/ I \ -2+1/21/ 1 



A 



^ fie f2eB£\ I 2eBt\ ^ , . . „^ e^fi f eB\ 



(33) 

where g{-) is a dimensionless function that is regular at 1. The analytic terms give a 
generic expansion, with the constant term representing, for instance. Landau diamagnetism. 
This piece includes contributions that have not been captured by the poles in ([t]), as this 
formula has zoomed in on the low energy states near the Fermi surface. The second term 
is the leading contribution coming from the Fermi surface and agrees with the previous 



computation (32). From (33) we can see that the oscillating term strongly dominates the 
susceptibility for | < z/ < ^. 

Finally, we can check that the scaling ([3]) is potentially observable in an experimentally 
interesting regime without being exponentially suppressed by temperature. Setting all 
dimensionless parameters except for v to be order unity, we can estimate the magnitude 
of the oscillations. Taking fi to be of order eV, T to be of order Kelvin and reinserting 
fundamental constants the exponent in our final result ([T]) is of order 

where Fb is the frequency of the oscillations measured in Tesla. In measurements on the 
underdoped cuprates, for instance, Fb/B ~ 10 [3], and so the exponent is not too large for 
a wide range of values of u. 



5 Discussion 

Using the holographic correspondence we have obtained the amplitude of quantum os- 
cillations in a family of strongly interacting quantum critical theories. Our expression ([T| 
provides a theoretical template for possible violation of Lifshitz-Kosevich scaling of the am- 
plitude with temperature due to strong interactions. We also found that at the marginal 
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value of the critical exponent ^ = ^, the Lifshitz-Kosevich result ([2| survives the inter- 
actions. Our results are perhaps the most concrete yet to emerge from applications of 
holography to condensed matter physics. The scalings we have described could conceivably 
be found in systems of current experimental interest. 

It will be important to generalise out computations to include disorder and to see to 
what extent the textbook Dingle scaling is modified. The dynamics of holographic theories 
with disorder has barely been studied [29]. Furthermore, while the singular loci ([T]) for 
the Green's function is the simplest following from the holographic correspondence [22], it 
is likely not unique. As new finite density dual geometries become available, it will be of 
interest to see to what extent our result (fll) is modified. 
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